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Abstract 



Oh 



Underground searches for dark matter involve a complicated interplay of particle physics, nuclear 
physics, atomic physics and astrophysics. We attempt to remove the uncertainties associated with 



astrophysics by developing the means to map the observed signal in one experiment directly into 

a predicted rate at another. We argue that it is possible to make experimental comparisons 

^ that are completely free of astrophysical uncertainties by focusing on integral quantities, such as 

> 

g(v m in) = f v . dvf(v)/v and f Vth h dvvg(v). Direct comparisons are possible when the v m i n 

Q\ space probed by different experiments overlap. As examples, we consider the possible dark matter 

signals at CoGeNT, DAMA and CRESST-Oxygen. We find that expected rate from CoGeNT in the 

XENON10 experiment is higher than observed, unless scintillation light output is low. Moreover, we 

determine that S2-only analyses are constraining, unless the charge yield Q y < 2.4electrons/keV. 

^ For DAMA to be consistent with XENON10, we find for qjy a = 0.3 that the modulation rate 

H 

must be extremely high ( > 70% for m x = 7GeV), while for higher quenching factors, it makes an 
explicit prediction (0.8 - 0.9 cpd/kg) for the modulation to be observed at CoGeNT. Finally, we find 
CDMS-Si, even with a 10 keV threshold, as well as XENON10, even with low scintillation, would 
have seen significant rates if the excess events at CRESST arise from elastic WIMP scattering, 
making it very unlikely to be the explanation of this anomaly. 

PACS numbers: 
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I. INTRODUCTION 



The search for dark matter is a central priority of modern high energy physics. Un- 
derstanding the nature of the ~24% of the universe composed of dark matter would give 
insight into the origin of galaxies and cosmic structures, the universe at high temperatures 
and a broader particle physics context for the Standard Model. To this end, a wide range of 
underground detectors have been and will be deployed to search for it (e.g., [IHS]). These 
programs utilize many different technologies and targets, and have sensitivities to a number 
of different WIMP scenarios, as well as mass ranges and interaction properties. 

These rare-event searches are some of the most sensitive detectors ever built, seeking a 
signal that may be as small as a few events per year. Consequently, they are sensitive to new 
and unexpected backgrounds as well. Any claim of dark matter discovery must be confirmed 
with multiple technologies before it can be believed. At the same time, comparing different 
experiments is a great challenge, with uncertainties from particle physics, nuclear physics, 
atomic physics and astrophysics compounding one another. 

Particle physics uncertainties can be explored by considering phenomenologically varied 
models, such as spin-dependent \±0\ or momentum-dependent [HJ [12] couplings, inelastic 
[13] [H] scattering, electromagnetic charge radius or dipole interactions [T5l - f2~T] . resonant 
dark matter [22], mirror matter [231 121] among others. Atomic and nuclear physics un- 
certainties can be better understood through theory and careful experiment (e.g., [2"oT - |2"T] ) . 
Astrophysical uncertainties are a challenge, however, for a number of reasons. 

High resolution numerical studies [28H30] have confirmed that Maxwellian distributions 
are generally good approximations to the phase space structure of dark matter halos. 
Nonetheless, significant deviations are found at high velocities that can be relevant for some 
scenarios, such as light WIMPs or inelastic WIMPs [291 SB E2], and in analytic solutions 
to NFW profiles such deviations are calculable [33]. Although the resolution of the numer- 
ical studies is not adequate to study small-scales accurately, variations from place to place 
within a halo can be important [29]. Pronounced structures such as streams or subhalos 
can dramatically alter expectations [MUST] . And of course, the "unknown unknowns" are 
impossible to quantify. 

It is important therefore to find methods that constrain scenarios without appealing to 
any model of the dark matter distribution. Some efforts at this have been studied already. 
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For instance, [38] argued that an independent comparison for the iodine spin-independent 
explanation of DAMA could be made by studying the comparable range of energy at a 
Xenon target, given their kinematical similarity. It was pointed out in that there is an 
overlap in velocity space between the ~ IkeVee signal at CoGeNT and the 7 keVr threshold 
at CDMS-Si. With positive results at two experiments, a measurement of the WIMP mass 
can be done without assuming a halo model [10]. Finally, [JI] studied the possibility of 
extracting f(v) from dark matter experiments in the future when large signals have been 
found. 

In this paper, we take a different approach. Rather than attempt to find the physical 
function f(v), or study variations in it, we attempt to directly map experimental signals from 
one detector to another. We do this by focusing on integral quantities, namely g{v m i n ) = 
( dvf(v)/v and f dvvgiv). We determine the robustness of constraints by considering 
the relationship between recoil energy and v min space, rather than actual velocity space. 
Although in our approaches we will gain less information about astrophysics, we can compare 
experiments even when f(v) cannot be reliably extracted. 



II. v min RANGES AND ASTROPHYSICS-INDEPENDENT SCATTERING 
RATES 

Our approach will be simple: we will endeavor to map an energy range in a given ex- 
periment into the halo velocity space, and from there into any other experiment we wish to 
compare to. In this way, we can determine what energy ranges of experiments can be di- 
rectly compared. In optimal situations, we will be able to extract g(v), while in less optimal 
situations we will only be able to discuss total rates. 

We begin with the differential rate at a direct detection experiment, which for elastically 
scattering DM is given by, 

dR N T M T p 

dW R = ^? a{ER)9{v ^^ (1) 

where \i is the DM-nucleus reduced mass, and Nt = kNatUp/Mt is the number of target 
scattering sites per kg with Na Avogadro's number and k the mass fraction of the detector 
that is scattering DM. The function g(v m i n ) is related to the integral of the DM speed 
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distribution 1 , f{v,t), by, 



oo 



g(Vmin,t) = I d,V ^ V,t ■ (2) 



There is a minimum speed that the DM must have in order to deposit recoil energy Er in 
the detector. For elastically scattering WIMPs this minimum velocity is 



'M T E R 

Making a comparison between different experiments is confused by the fact that it is 
not a single velocity that contributes to the scattering rate at a particular Er. Rather, all 
particles with velocities greater than v min will contribute, making it impossible to map rates 
into velocity space. 

However, we can consider a related space - t> miri -space, whose elements are the sets of 
all particles with velocities greater than v m i n - Because all particles with adequately high 
velocities contribute, it is reasonable to consider a mapping between Er and v m i n through 

©■ 

This simple relationship allows us to compare results from different direct detection ex- 
periments without making an assumption about the distribution of DM velocities in the 
Milky Way's halo, provided one can relate the scattering cross sections at the various ex- 
periments. In the standard cases of SI or SD DM the nuclear scattering cross section can 
be related to the nucleonic (in this case the proton) cross section as 



/i 2 (f p Z + f n ( A-Z)) 2 
>p u 2 f 2 



asi(E R ) = a p ^- " JP ~ " F {Er) (4) 



>2 



*so{Er) = ^y^t , (5) 



a p /i 2 {a 2 p S pp {E R ) + a p a n S pn {E R ) + a 2 n S nn {E R )) 

allowing comparison of different experiments, we have defined /x nx as the DM-nucleon re- 
duced mass. 

Let us suppose we have two experiments to compare, with targets Ti,2 with masses M\^- 
We assume the first has a signal which appears over an energy range [E^, E^ h ]. This 
energy range correspond to v min ranges [v l ^ n ,v^}, using M. 



1 It is usually assumed that the DM follows a Maxwell-Boltzmann distribution (in the galactic frame), with 
characteristic speed vq, in which case (again in the galactic frame) f(v) oc v 2 e~ v l v ° . 
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This brings us to the central point of our efforts: to make a comparison between two 
experiments one must first determine whether the v m i n space probed by the two experiments 
overlaps. As a matter of practical course, a given experiment has a lower energy threshold 
E m in, which can be translated into a lower bound on the v m i n range. If experiment 1 has 
data for the differential rate of DM scattering in their experiment, dR\jdE R at energies 

(2) 

this can be used to predict a rate at energy E\ at experiment 2, dR 2 /dEji, or vice versa if 
experiment 2 has the signal. Thus, we have 

[■^low) E high] ^ [ V mini V min\ ^ [E\oL > ^hiq h\ > (6) 



where 



W (2) E (2) , _ f4M™ (1) (1 ) 

L-^Zoio' 111 high! ~ „ , ,(21 L^oiu) ^highV \ ' ) 



We can invert ([!]) to solve for g{y min ) limited to the range v min E [v. 



low high I 
min,l' ^min,l\ 



2m x /i 2 (i-Ri , . 

JMAK-rrip pa{hji) dhi 



This then allows us to explicitly state the expected rate for experiment two, again 2 re- 

?( 2 ) 1 
highi ' 



stricted to the energy range dictated by the appropriate velocity range i.e. EE [El o ^E h 



Analogous to the energy mapping above, we have a rate mapping, 



with 



dR\ , \ dR 2 

dE, ^ 9{Vmm) ^ dE^ (9) 



dR 2 _ a 2 {E 2 ) dR t ( ix\Mf 

dE R [ 2) (#M? P \dE R [fi M p 21 



Equations 0, ® and goj are the central results of this paper. They make no astrophysical 
assumptions, but only rely upon the assumption that an actual signal has been observed. 

We now focus on the SI case, since there are a greater number of experiments probing 
this scenario, but the analysis for SD is similar. In this (SI) case we can use (Q to rewrite 



2 Since g(v), by its definition, is a monotonically decreasing function of v m i n , one can in principle go to 
lower energies as well, but one may only place a lower bound on the predicted rate, rather than make a 
true prediction. 
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(10) in a simple form 



dR 2(F ,_ 4 2) W dRi ( \A M? \ (u , 
dEn ™ <#> 2 \ dE R ^1M« V ' 1 J 

1 2 J 

where we have introduced a target specific coefficient 

C^ = n^(f p Z^ + f n (A^-Z^)) 2 . (12) 

In certain situations differential rates may not be available and instead it is only possible 
to compare total rates, this is the situation at present with CRESST. In general the total 
rate at a particular experiment with energy — and corresponding velocity — thresholds of 

(Ei ow ,v l ™ n ) and (E high ,v^), can be expressed as, 

2N Anm k /""high 
R = Z -^P^^- / dve(E R )a(E R (v))vg(v) . (13) 

m X M T Jv low 

For the particular case of SI on which we are focused this becomes, 

where e(E R ) an an energy-dependent efficiency. To compare two experiments, we must 
extract the energy dependent terms from the integral. So while we make no assumptions 
about g(v), we evaluate the form factor at a value E 2 = Ei/i^M^' / where the ra- 

tio £2(^2) E2(E 2 )/ei(Ei)E 2 (Ei) is minimized or maximized, depending on whether we are 
considering a putative signal or constraint. Thus comparisons of rates at two experiments 
may then be simply compared by taking ratios of Ct with the form factor evaluated at the 
conservative value E, 

e 2 (E 2 )Fj(E 2 ) 

In order to determine what comparisons can be made between experiments, we must ex- 
amine the relevant velocity space they probe. We re-emphasize that the signal at energy 
E[ ow < E < Ehigh is sensitive to all particles with velocity greater than v m i n (E, M^, M x ) 
through the integral g(v m i n ). A separate experiment with threshold E will offer constraints 
independent of astrophysics if the resulting minimum velocity v < v^fn i- The optimal limits 
are reached when v < v l ™ nl . We illustrate this in Fig. [I] for an ensemble of experiments, 
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FIG. 1: v m in thresholds for various experiments. Solid bands are CRESST Oxygen band, 15- 
40 keV (red, top), DAMA Na band 6.7-13.3 keV (green, middle), CoGeNT Ge 1.9-3.9 keV (blue, 
bottom). Constraints are Xenon 1, 2 and 5 keV (dashed, dotted, and dot-dashed, thick blue), and 
CDMS-Si 7 and 10 keV, (dot-dashed and dashed, thin red). 

some with signals, some without. The possible comparisons between these various exper- 



iments will be the subject of the subsequent sections. Using (11) scattering rates can be 
compared between experiments. However, to compare to actual experimental data the rel- 
ative exposures, efficiencies and other detector-specific factors must be correctly taken into 
account. In the next section we describe in detail the experimental parameters necessary 
for the comparisons in the rest of the paper. 

III. APPLICATIONS: A COMPARISON OF EXISTING EXPERIMENTS 



The important consequences of (10) are immediately obvious. In principle, one can com- 
pare a positive signal at one experiment with one at another, or test the compatibility of a 
null result with a positive one. Unfortunately, ideal circumstances will rarely present them- 
selves: additional backgrounds can complicate the extraction of g(v), resolution can smear 
signals, or uncertainties in atomic physics (such as quenching factors) can complicate issues, 
making a precise extraction of the true E^r and hence v m i n impossible. Furthermore, the 
signal may appear as a modulation (as in DAMA) limiting access to g(v) to a summer/winter 
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difference, or a lower bound on its mean. Finally the signal may be of such low statistics 
that a reliable inference on the shape of g(v) will be impossible, as is expected in many 
experiments before scaling to larger targets or running for longer exposures. 

Nonetheless, in light of these challenging issues, there remain meaningful comparisons 
that can be made between experiments. Especially since these transformations preserve all 
information in the signal, we should be able to make the strongest possible relative state- 
ments without invoking additional assumptions about the halo. Such results are especially 
interesting in view of recent results that may pertain to light WIMPs. Since light WIMPs 
probe the highest part of the velocity distribution, where deviations from Maxwellian prop- 
erties are the most likely, our approach is especially relevant. 

We consider three potential signals: the CoGeNT low-energy excess [TJ , the DAMA annual 
modulation signal [42] and the recently reported Oxygen-band events at CRESST [43J and 
for constraints: XENON10 (both conventional analyses and S2-only) and CDMS-Si; we 
describe the relevant parameters necessary for comparison between the various experiments 
below. 

CoGeNT 

The CoGeNT experiment [TJ consists of a low noise germanium detector with 330 g of 
fiducial mass which has reported data for 56 days of exposure. CoGeNT reports recoil 
energies that range from ~ 0.4 keVee to ~ 12 keVee, but we focus here on the events 
between 0.4 keVee and 3.2 keVee. The observed electron equivalent energy is related to 
the nuclear recoil energy by, E obs = 0.2(i? r /keV) 112 , so that the range of nuclear recoils 
of interest is 1.9 — 12keV. In this range there are two cosmogenic peaks whose position 
and width are well understood, a relatively flat spectrum above these peaks and a clear 
excess at energies below the peaks. It is this low energy excess that may be due to a DM 
signal and, rather than assume a particular functional form and fit, we extract it from 
the data by taking the data below the first peak (E ee < 1 keVee) subtracting from it the 
average of the high bins (E ee > 1.6 keVee). Thus our "signal" region, shown in Fig. [Tj is 
0.42 keVee (1.9 keV) < E ee (E r ) < 0.92 keVee (3.9 keV). Finally, we take into account the 
detector efficiency [TJ. 
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DAMA 



DAM A, which has a Nal target, has accumulated a total of 1.17 ton-years of data, from 
both DAMA/Nal 02] and DAMA/LIBRA HUE]. They have observed, at the 8.9 a C.L., an 
annually modulating signal, in the 2-6 keVee energy range, whose phase is consistent with 
that expected from DM. For low DM mass, recoils off the sodium dominate the spectrum, 
we take the quenching factor to be qNa = 0.3, although we will discuss the effects of varying 
this. We concentrate here on the low energy range of recoil energies, 2-4 keVee, and in Fig. [T] 
we show the f min -space for this range of energies, assuming a quenching factor q Na = 0.3. 
A lower (higher) q^ a will result in the band moving higher (lower) in the plot. Specifically, 
for qN a pa 0.45 CoGeNT and DAMA are probing nearly identical ranges of velocity space. 

In addition to the modulated rate DAMA has measured the total rate of recoil events and, 
as emphasized in [38], the DAMA unmodulated rate may also provide non-trivial constraints 
on models of DM. We do not consider the effects of channelling, which are believed to be 
small [IS]. We take into account the energy resolution at DAMA, by smearing with a 
Gaussian distribution of width <x(E)/E = 0.448/v 7 ^ + 0.0091, with E in keVee. 

CRESST 

CRESST consists of 9 CaW0 4 crystals (and 1 ZnW0 4 crystal). They recently reported 
[43] an excess of O-band events in approximately 400 kg days of exposure. Although they 
have not yet reported a detailed spectrum for each crystal it is still possible to gain some 
spectral information from Ref. [43j. 

The upper bound on the energy of the events is set by their 40 keV upper limit on their 
search box. For the lower bound we use a threshold of 15 keV. In reality, the threshold in 
each of the nine detectors is determined by the value where the leakage is expected to be 0.1 
events, which ranges from 9.65 to 22.65 keV. In this total range 32 events are seen with an 
expected background of 8.7. In Ref. [43] individual detectors are listed and two detectors are 
explicitly plotted, and consequently we can determine that for seven detectors, there are 22 
of 27 events above 15 keV, with an expected background (in these detectors) of 7.2 events, 
where to be conservative we have attributed all neutron and gamma backgrounds to these 
seven detectors. Thus, since we only use data from 7 detectors, and taking into account 
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an efficiency of 90% (as used in the commissioning run), the data considered here has total 
exposure 280 kg days, and 22 events between 15-40 keV, with an expected background of 
7.2. In the future, when all thresholds, exposures, and events are reported by the CRESST 
collaboration, these results can be refined, but for the moment, we make these conservative 
assumptions. The range of threshold velocities (i.e., v m i n values) corresponding to CRESST 
O-band events between 15-40 keV is shown in Fig. [1} 

For constraints, we consider those experiments that are particularly sensitive to light 
WIMPs: XENON10 (both a conventional analysis with a range of scintillation efficien- 
cies, as we shall describe, as well as S2 only) and CDMS-Si. 

XENON10 

We use the unblind XENON10 analysis jlH] on a 5.4 kg active target of Xenon taken over 
58.6 days between October 2006 and February 2007. This analysis found 13 events between 
16 keV and their upper threshold of 75 keV. Their lower threshold is set by requiring a 
minimum of ~ 12 ionization electrons in the S2 signal, which for a constant £ e // = 0.19 
corresponds to a threshold of ~ 2keV, in addition there is an analysis threshold on the SI 
signal of ~ 5 keV, with the same assumption on £ e ff- We will also consider a potential 
S2-only low-threshold (7 drift electron ~ 1 keV) experiment as recently discussed by [47], we 
consider the efficiency adjusted exposure to be 5.1 kg days. For the charge yield considered 
in [47], this would correspond to a threshold of 1 keV. In Figure [I] we show limits for a 
Xenon experiment (thick blue lines) with a 5 keV threshold (upper line), as well as a what 
sensitivities 1 and 2 keV (lowest line and middle line, respective) thresholds would achieve. 

There is considerable uncertainty in the behavior of £ e //, especially at low energies, and 
here we will consider the three cases discussed in [IS] which we denote by MIN, MED and 
MAX ordered by increasing value of £ e // at recoil energies of 2 keV. The corresponding 
detector resolutions and efficiencies are taken from |48j . 

CDMS-Si 

The CDMS experiment contains both germanium and silicon detectors, we focus here on 
Si since it is sensitive to lighter WIMPs due to smaller mass and lower thresholds. There 
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have been several analyses of silicon data taken at the Soudan mine [4*9l451~] which combined 
have a raw exposure of 88.6 kg days. For these analyses we use the efficiency presented 
in [5T] which has a threshold for nuclear recoils of 7 keV. In all data there are no signal 
events observed below 50 keV. We do not consider here the recent low threshold analysis of 
CDMS-SUF [52]. 

In Fig. [T] we show the regions of velocity space associated with the experiments discussed 
above, the bands denote the range of v m i n necessary to see events in the experiments with 
possible DM signals and the region between the curves are the same quantity at various null 
experiments. We show both 7 and 10 keV thresholds for a silicon experiment, we choose 
these limits as at 7 keV the CDMS efficiency is going to zero, making limits that are reliant 
on that threshold more sensitive to uncertainties, while at 10 keV it is stable at ~ 20%. For 
a Xenon target we consider 1, 2 and 5 keV thresholds, which may be reachable for different 
assumptions regarding C e ff and Q y . In Table [I] we show, for fixed m x = 10 GeV, how the 
energy ranges probed at CoGeNT, DAMA and CRESST translate to various elements used 
in other experiments. 

A careful examination of Fig. [T] shows a number of things: first, CoGeNT can be tested 
only by the 7 keV threshold of CDMS-Si, as well as a Xenon analysis sensitive to low ener- 
gies. The first two have sensitivity below 10 keV, while Xenon can only make astrophysics- 
independent statements if the threshold is lower than 2.5 keV. This demonstrates, explicitly, 
that a model-independent comparison involves reaching signals present at 2 keV, and clar- 
ifying the scale to which £ e // must be measured. The CDMS-Si analysis is only sensitive 
right at its threshold to CoGeNT. As questions have been raised [53] about the precise value 
of the Si threshold, if one restricts oneself to the higher threshold, no limits are possible. 

Similarly, DAMA is tested well by CDMS-Si with a 7 keV threshold, but only marginally 
at 10 keV. Low-threshold Xenon analyses can give robust limits of DAMA, while higher 
thresholds are generally limited to heavier masses. 

Finally, we see that the CRESST results are completely tested by the low-threshold 
XENON10 analysis, CDMS-Si (even with a 10 keV) threshold. While the nominal threshold, 
depending on the details of £ e// , of XENON10 (~ 5 keV) and XENON100 (~ 6 keV) is too 
high, both experiments can probe down to 4 keV with moderately reduced sensitivity, and 
energy smearing will given XENON sensitivity to the CRESST signal. 
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Approx. range 


Na 


Si Ar 


Ge Xe 


CoGeNT (Ge): 2 


- 4 


4.3 - 8.6 3.9 - 7.8 


3.6 - 7.2 3.0 - 6.0 


2-4 1.3-2.5 


DAMA (Na): 6 - 


13 


6.6 - 14 6-13 


5.5 - 12 4.6 - 10 


3.1 - 6.7 1.9 - 4.2 


CRESST (0): 15 - 


40 


15 - 40 14 - 36 


12 - 33 10 - 28 


6.9 - 19 4.3 - 12 



TABLE I: Conversion of energy ranges (all in keV) between various experiments/targets for a 10 
GeV DM particle, using the expression in Q. 

With these ranges in hand, we can proceed to compare the experiments directly. We 
shall see that if the potential signal is large enough, g(v) can be extracted directly, even if 
f(v) cannot be extracted with any reliability. In such cases, we can make slightly stronger 
statements involving the spectra. However, even if g(v) cannot be reconstructed, we can 
still make significant statements by integrating over the relevant velocity range. 



A. Application I: Employing Spectra in Near-Ideal Situations (CoGeNT) 

We consider first the situation when there is sufficient data to be able to extract a recoil 
spectrum, CoGeNT is a example of such an experiment, because the putative signal is quite 
large. We concentrate on the events below 3.2 keVee where the DM signal should be largest 
and there are few cosmogenic backgrounds. In this range, in addition to the possible DM 
signal at low energies, the data contains several clear cosmogenic peaks and a constant 
background above the peaks. We average the [1.62-3.16 keVee] bins as an estimate of the 
constant background and subtract this from the bins in the [0.42-0.92 keVee] range, which 
we then consider as the DM signal, after this subtraction there are 92 signal events before 
efficiency correction. This allows us to determine g(v) or, equivalently, predict the rate 
at any other experiment in the equivalent energy range. One can easily observe from its 
definition that g(v) is monotonically decreasing as a function of v (see, for instance the 
discussion in |4"I]). and thus the value at the low end of this range is a lower bound for lower 
values of v. This is not especially relevant for our analysis here, but would be likely relevant 
in situations where the other experiments could probe lower energies as well. 

Since we will compare this with the XENON10 experiment, we choose f p = 1 and f n = 0, 
which is motivated from light mediators mixing with the photon, since it will give the most 



lenient bounds. Using (11) we can map the CoGeNT signal onto a Xenon target, and study 
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FIG. 2: The extracted CoGeNT signal (left and bottom axes) and the rate it is mapped to on a 
Xenon target (top and right axes) for m x = 10 GeV (rescaled by form factors at the corresponding 
energies F^- e (E^ e ) , FQ e (E^ e ) ~ 1). The dashed line is the lower bound on the rate at low energies, 
using the monotonically falling nature of g(v m i n ). 

the signal that would arise at XENON10. We show this in figure [2] 

What is remarkable about this figure is that - once the CoGeNT signal is specified - the 
expected rate on a Xenon target is completely unambiguous (and similarly on any other 
target). This involves no assumptions about the halo escape velocity, velocity dispersion, or 
even the assumption that the velocity distribution is Maxwellian, but requires only an input 
of the WIMP mass. 

After taking into account exposure and the detector efficiencies (MIN, MED and MAX 
cases described above) we can predict the total number of events predicted by the CoGeNT 
events (if they are indeed coming from elastically scattering DM), we show this in Fig. [3j 
Since there were no events at XENON10 in the energy range corresponding to the CoGeNT 
range we see that independent of all astrophysical assumptions, only for L^jj N are CoGeNT 



o C.L. In the MIN case, m x < 11 GeV allows Co- 



and XENON10 are consistent at the 
GeNT to evade XENON10. For MED and MAX cases the predicted signal at XENON10 
would be too large by a significant amount, excluding the elastic SI WIMP scattering inter- 
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FIG. 3: The number of events predicted at XENON10 by the possible DM signal at CoGeNT for 
3 cases of C eff , MIN (dashed red), MED (solid green) and MAX (dotted blue). The black line is 
the 90% C.L. upper limit on the number of events allowed by XENON10 data, the region above 
this line is excluded at 90% confidence. 

pretation by more than an order of magnitude. 

Because of the uncertainties associated with extraction of the value of C e ff at low energies, 
additional attempts have been made to probe the low energy region with Xenon experiments. 
In particular, [171 [M] examined data from XENON10, and used only the ionization signal 
(S2), which is typically larger than SI and can allow a more reliable signal at low energies. 
The value of the charge yield (drift electrons per keV) was extracted from Monte Carlo. 
Using the values there, the equivalent energy range for CoGeNT is approximately 8 ~ 13 
electrons, above the 7 electron threshold. Assuming a value of Q y = 4 electrons/keV for 
instance, the threshold of 7 electrons at XENON10 only captures a portion of the signal 
predicted by CoGeNT. 

While the 7 electron cutoff corresponds to a particular value of energy in principle, Poisson 
fluctuations smear this. Nonetheless, an interesting question is the expected rate on the 
target used by (UJ |5l], with 5.1 kg d of effective exposure. This is most easily phrased in 
terms of the question of what charge yield can make these experiments consistent. Assuming 
a constant charge yield over the energies in question, we can calculate the likelihood based 
on Poisson fluctuations of events appearing in the XENON10 experiment, which we show in 
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FIG. 4: (left) The number of events predicted (labels on contours), by CoGeNT, at XENON10 for an 
S2 only analysis |47} 154j for various S2 thresholds, assuming a constant value Q y = 4electrons/keV. 
(right) The signal above threshold of 7 electrons, but assuming different constant values of charge 
yield, Q y . 

Figure [ij One sees that one would require a charge yield of roughly Q y < 2.4electrons/keV 
for consistency, much lower than the value of Q y ~ 7 extracted by [HJ Whether such 
a significant difference is reasonable will no doubt be subject to a great deal of discussion 



B. Application II: Total Rate Comparisons in Sub-Optimal Situations (CRESST) 

The above situation with CoGeNT is close to ideal: low backgrounds, high statistics, good 
energy resolution and calibration. In contrast, there are often situations with significantly 
less ideal characteristics. In particular, it may be that not enough is known about the 
backgrounds, or the data itself, to be able to extract a recoil spectrum for DM, but we shall 
see it is nonetheless possible to say something about the total number of DM scatters. This is 
the case for the CRESST data, which we estimate has 15 events above background between 



15 and 40 keV (see the discussion in III). We use (15) to compare the CRESST integrated 
rate to the null results of both CDMS-Si and XENON10, Fig. [5} When comparing the two 
experiments we take into account efficiencies and form factors so as to be as conservative as 
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FIG. 5: LH plot: the CRESST prediction for the total number of events at CDMS-Si (solid red) 
and XENON10, for C^/ f N (dashed black) and £^ f f D (dotted black), the dotted (horizontal blue) 
line is the 90% C.L. upper limit on the number of events allowed by CDMS-Si, the region above 
this line is excluded at 90% confidence. RH plot: the 90% C.L. upper limit on the number of 
events at CRESST as predicted by CDMS-Si (solid red) and XENON10, again for Cf/ f N (dashed 
black) and £-eff D (dotted black), the dotted (blue) line is the number of events we estimate above 
background in CRESST. 



possible, as explained after (15) 



As is clear from Fig. [5] any sizeable signal in this range is highly incompatible with both 
the XENON10 and CDMS-Si results. While some have criticized the calibration at the 
lowest energies for CDMS-Si [S3], the lowest energy relevant for 15 keV Oxygen recoils is 
above 10 and typically 11 keV on Silicon, depending on the WIMP mass. Thus, these 
constraints are likely quite stable to future modifications, making elastic WIMP scattering 
very unlikely to be the explanation of the CRESST anomalous events. 



IV. OTHER APPLICATIONS AND FUTURE RESULTS 

DAMA also has extracted a recoil spectrum, possibly associated with DM, but in this 
case it is for the modulating part of the DM signal, i.e. DAMA allows extraction of g(v,t). 
We can repeat the exercise of translating from one experiment to another to get a prediction 
for the size of the modulating signal at XENON10. Since XENON10 took its data in the 
winter and saw no events in the region corresponding to DAMA's 2-6 keVee, this places an 
upper limit of 2.3 events in the winter which in turn places a lower bound on the amount of 
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FIG. 6: The 90% C.L. lower limit on the modulation fraction allowed by XENON10 data, for a 
quench factor in sodium of 0.3 (LH plot) and 0.45 (RH plot) and for 3 cases of C e ff, MIN (dashed 
red), MED (solid green) and MAX (dotted blue). A modulation fraction below the curves is ruled 
out at 90% confidence. 



modulation the DM signal must have in order not be ruled out by XENONlO's null result. 
We present this lower bound on the modulation fraction 3 in Fig. [6] for two choices of the 
quench factor in sodium, q^ a = 0.3, 0.45. Thus, irrespective of astrophysics, in order for 
DAMA to be consistent with XENON10 the modulation fraction has to be larger than 20% 
and in most cases almost 100% for the standard assumption of gv a = 0.3. For the more 
extreme choice of q^ a = 0.45 the modulation may be smaller but for DM heavier than 10 
GeV it again has to be above 20%. 

An interesting relationship between CoGeNT and DAMA can be made here. The modula- 



tion at DAMA can be applied to CoGeNT through (15). In doing so, one finds a modulation 



O(0.8 — 0.9cpd/kg) expected at CoGeNT. With a quenching factor q^ a = 0.3, this is ex- 
pected to overlap the L-shell peaks, which, in decaying away, would make a rising signal 
difficult to extract. The modulation in the signal range we cannot predict. 

On the other hand, if q^ a — 0.45, then the energy range of CoGeNT overlaps that of 
DAMA. The 0.8— 0.9cpd/kg modulation amplitude would then be visible over the ~ 5cpd/kg 
in that range (i.e., a modulated amplitude of ~ 20 %, or ~ 40% peak-to-peak), which should 
be detectable over an annual cycle. 

Finally, it is intriguing to employ this technique to study what sensitivities future exper- 



We define the modulation fraction as | + ^ where S, W denote the summer and winter event rate respec- 
tively. 
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iments will have to existing signals. In particular, we can consider the COUPP experiment, 
with a CF 3 I target. Focusing on scattering off of Fluorine, the CoGeNT events would be 
visible in a range 4.7 - 9.4 keVr, while the DAMA modulation would be present in the range 
7.3 - 14.6 (4.8 - 9.7) keVr for q Na = 0.3(0.45), for m x = 7 GeV. Thus, a threshold in the 5 
- 9 keVr range should allow astrophysics independent tests of these signals. 



V. DISCUSSION 



The search for WIMP dark matter has as its elements three central goals: to discover 
the WIMP, and to measure its mass and interaction cross section with matter. Although we 
have proceeded for years without a confirmed discovery, the focus has been on what ranges 
of mass and cross section are excluded. Regrettably, this thinking has crept into our whole 
approach to discussing a comparison of WIMP searches — we compare compatibility within 
the confines of m x — a plots, confusing the answer to the latter two questions (the WIMP 
mass and properties) with the central first question: has dark matter been discovered? 
To do so necessarily entangles our astrophysical uncertainties into our results, and even 
worse, makes it difficult, if not impossible, to determine how sensitive the conclusions are 
to variations in the halo model. 

To this end, we have explored a new technique to compare experiments in the presence 
of a positive signal, and importantly, to do so without invoking any astrophysical model 
whatsoever. The energy range of a given experiment can, for a given WIMP mass, be 
mapped unambiguously into an equivalent energy range at another using the expression in 
(J7]), allowing an apples-to-apples comparison of signals and limits. Moreover, the measured 
rate at one experiment can be mapped into a rate at the other experiment with the expression 



in (10), once the particle physics model is specified, yielding a completely unambiguous 
prediction for the second experiment. 

This is done by implicitly solving for the function g(v). In a sense, each experiment 
is actually a measurement or upper limit of the function g(v). This motivates a new and 
simple comparison of experimental results by simply showing the different values and limits 
extracted for this function from different results, which we do in Fig. [7} 

To determine this plot, in the presence of a positive signal, one needs merely to read off 
g(v) from ([8]). In the absence of a (clear) signal, there is always a certain element of choice in 
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FIG. 7: A comparison of measurements and constraints of the astrophysical observable g(v) [see 
relevant expressions in 0,0) ,@] for m x = 10 GeV: CoGeNT (blue), CDMS-Si (red, solid), 
CDMS-Ge (green, dot-dashed), XENON10 - MIN C eff (purple, dashed), and XENON10 - MED 
C e ff (gray, dotted). CoGeNT values assume the events arise from elastically scattering dark 
matter, while for other experiments, regions above and to the right of the lines are excluded at 
90% confidence. The jagged features of the CDMS-Ge curve arise from the presence of the two 
detected events. 

how one quantifies a constraint. However, one can exploit the fact that g is a monotonically 
decreasing function, so for our constraints, we simply assume that g(v) is constant below 
v, and assume a Poisson limit on the integral of ^ from the experimental threshold to v. 
However, other techniques could also be used, see the Appendix for more details. 

This approach with a g — v plot has numerous advantages over the traditional m x — a 
plots. It makes manifest what the relationships between the different experiments are in 
terms of what t> m j n -space is probed, and shows (for a given mass) whether tensions exist. 
Moreover, the quantity g(v) is extremely tightly linked to the data, with only a rescaling 
by form factor as in ([8]). Thus, unlike m x — a plots, which have a tremendous amount of 
processing in them, this provides a direct comparison of experimental results on the same 
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plot. 

In light of this, we would propose that all future results take a three-tiered approach to 
data comparison and presentation 

• m x — a plots: the current standard presentation is extremely useful for understand- 
ing implications for models irrespective of other experiments, and for making model- 
dependent comparisons of different results. We emphasize that our new techniques 
are complementary to this, and not a replacement for it. 



A mapping of contested or compared results: by employing (10) an experiment can 
state in a completely astrophysics-independent fashion whether two results are com- 
patible or not. The use of these tools would remove any discussion as to whether 
different assumptions about escape velocities would impact the question of consistency 
of conclusions. 



• A measurement of g(v) : By presenting results in a g — v plot as in Fig. [7j the astrophys- 
ical dependences of results become more obvious, and relative (in) consistencies clear. 
For experiments with only a single nuclear target, the plot can be made for a single 
value of m x , and can be unambiguously mapped to any other value. For experiments 
with multiple targets, this will be slightly more involved, but still will likely require 
plotting the result for only a few values of m x . 

In the asymptotic future, we certainly hope for a large number of experiments, with 
different targets and technologies, all seeing results with high statistics. In the meantime, 
the question of whether dark matter has been discovered prompts the need for new techniques 
that do not rely on outside assumptions, but at the same time, keep the maximal amount 
of information possible. We have proposed such a technique, by mapping signals from one 
experiment to another via (10), comparable energy ranges through ([7]), and by extracting 



the physical astrophysical quantity g(v) through ([8]). 

We have shown numerous applications of this to existing results from CoGeNT, DAMA 
and CRESST. In particular, we have found that only for low values of £ e ff, a positive 
signal at CoGeNT can coexist with null results from the conventional XENON10 analyses. 
For S2-only analyses, however, we have seen that consistency requires a charge yield of 
Q y < 2.4 electrons/keV, well below most current assumptions. By comparing DAMA and 



20 



XENON10, for q^ a = 0.3, we have shown the modulation rate is typically higher than 80%, 
which would be difficult to achieve for elastically scattering WIMPs over a large range of 
2 — 4 keVee, even with non-standard halo models. On the other hand, for q^ a = 0.45, 
the modulation can be lower, ~ 10%. However, in this case, the DAMA modulation maps 
into the CoGeNT energy range, predicting a clear modulation amplitude of ~ 20% in the 
low-energy range. Finally, we have seen how even in situations with large backgrounds and 
unclear spectra, as is the case of the excess events reported by the CRESST collaboration, 
constraints by CDMS-Si and XENON10 unambiguously test this, making the elastic WIMP 
scattering explanation of them very unlikely. 

While it remains to be seen which if any of these signals will turn out to be genuine 
signals of dark matter, these tools will provide a means to remove one of the important 
uncertainties in their comparison. Perhaps, with higher statistics, this can be inverted, and 
by requiring the consistency from different measurements of g{v), dark matter will yield 
information on astrophysics, as well. 
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Appendix A: Displaying direct detection results in an astrophysics independent 
fashion 

As well as allowing comparison of positive results between different experiments in a 
fashion that is independent of astrophysics, our technique also allows constraints to be 
compared to each other and to putative discoveries. Such comparisons can be made, such 
as we have done by constraining g(v) in Figure [7j Here we outline, in more detail, how this 
is carried out. 
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For positive results the comparison can be made at the spectrum level through the appli- 



cation of (10). This is most easily done in the situation that the statistics are large enough, 



and backgrounds low enough, that a meaningful rate dR/dE R can be extracted. Then, using 



(10), a direct measurement of g(y) can be given. In situations where rates are too low to 
simply read off g(v), alternative techniques would be needed. The simplest would just be 
to take a large enough bin in v such that statistics are adequate, but more sophisticated 
approaches, utilizing the monotonicity of g{v) would also be possible. We leave such studies 
for future work. 

If an experiment does not see a sufficient number of signal events to claim discovery, then 
it is likely that one will wish instead to place a constraint on the properties of dark matter. 
In general, one should first ascertain the bound on the parameterized WIMP cross section at 
the confidence level required. This can be simply done, using whatever confidence estimator 
is already used for astrophysics-dependent a — m x plots. 

Suppose that one wishes to employ some confidence estimator C(dR/dE R (m x )) to place 
limits, where dR/ dE R {m x ) is the expected recoil spectrum for some a , i.e., (JT]). This 
estimator may simply be using Poisson statistics, evaluating the integral of the spectrum, 
or using more advanced techniques that use spectral information as well, such as those 
of Yellin [56]. For a given value of m x , for instance, one varies other parameters until 
one achieves, e.g., C — 0.1 allowing one to claim a 90% exclusion for those parameters. 
Assuming that such an analysis has already been performed for explicit halo models, it is 
straightforward to place a bound on pcrg(v m i n )/m x , for a particular choice of DM mass, in 
the general astrophysics case. 

For standard a — m x plots, g(v) is fixed, for instance a Maxwellian distribution, with a 
fixed vq and v esc . The only free parameters in dR/dE R (m x ) are then m x and do (in the SI 
case), or a p and a n (in the SD case). In our case, since we do not want to use a Maxwellian 
g(v), we have an additional free parameter. 

Since g(v ) is a monotonically decreasing function an upper bound on its value at some 
velocity v\, g{v\) < g±, also applies to all lower velocities. Thus, the most conservative form 
that the upper bound on g{v) can take is that of a step function 

g(v;v 1 )=g 1 Q(v 1 -v) . (Al) 

Physically, this would correspond to stream in f{v) with velocity v\. 
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Using this, Q b 



ecomes 



dR 



N T M T p 



a(E R ) c/10 (ui - v min (E R )) . 



(A2) 



dE R 



2m x p 2 



For a given WIMP mass m x , the overall scaling is now proportional to e.g., pag\jm x in 
the SI case, rather than simply pa/m x as in the standard case where g is specified. For a 
given Vi, one can then place a limit on this combination using the existing estimator. 

In short: to calculate the appropriate limits on g(v), one should use whatever technique 
one was intending to use for the standard analysis, but now replace the Maxwellian g(v) 
with the step function form. For any given m x , one places a limit on pag\jm x as one would 
have on pa/m x , or, o for fixed p and m x , precisely as before. 
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